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ABSTRACT 

Analyses and interpretation of surtace layer and synoptic-scale data 
obtained in the Northeast Atlantic were performed to obtain descriptions 
of the evaporation duct and associated atmospheric and oceanic synoptic 
features. The surface layer data were quite unique because they were ob- 
tained from high quality measurements from ships spatially separated in a 
fixed array. Magnitudes and horizontal homogeneity of duct heights were 
compared to air-mass trajectories and weather patterns. The mean duct 
height was 4-5m and higher values (8-15m) occurred with North and North- 
west trajectories. During times when the heights were 8-13m, horizontal 
homogeneity occurred 35% of the time with a maximum duration of 12 hours. 
This was established on the basis of a two hundred kilometer separation 
between locations of duct height estimates. 

Significant features of the duct height and its temporal and spatial 
varlations were related to synoptic scale descriptions. Satellite imagery 
used in conjunction with point observations appears to provide the most 
useful information in describing the intensity and areal distribution of 


the evaporation duct. 
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I. INTRODUCTION 


Over the surface of the ocean there exists a phenomenon described 
as the evaporation duct. The military significance of its existence 
has been known for many years. Attempts have been made to quantify 
the evaporation duct height and determine its effect on electro- 
magnetic propagation. Assessment of evaporation duct characteristics 
in tactical military applications should only require routine ship- 
board measurements and simple calculations; hence, a great deal of 
study has been performed to parameterize the duct on the basis of 
Simple observational measurements of surface and near-surface para- 
Meets. 2150, a5 the use of low-flying supersonic aircraft and 
missiles become increasing threats to Naval surface combatants, timely 
and precise descriptions of the nature of the evaporation duct will be 
required to permit its exploitation. 

It has been expected that in an area in the high latitudes, well 
away from large land-mass influences, and where frequent migratory 
weather systems disturb the nature of the atmosphere near the surface, 
the height of the evaporation duct is likely to be small and its 
occurrence highly localized. Past research has been primarily limited 
to regions where ducts are known to be both strong (high) and per- 
Sistent (e.g. tropics, southeast portion of subtropical high, coastal 
regions of warm offshore winds) with relatively little examination of 
higher latitude, open oceanic regions. This study will examine the 


occurrence and nature of the evaporation duct over the North Atlantic 





Ocean subject to highly varying atmospheric conditions. This 
examination will focus on determining significant ducts and their 
1) frequency of occurrence and duration, 2) temporal and spatial 
nature, and 3) sensitivity to atmospheric conditions. 

The results reported here were made utilizing a unique set of data 
acquired from the Joint Air-Sea Interaction Project (JASIN) conducted 
in the Northeast Atlantic during the late summer of 1978. The results 
were analyzed and comparisons were made between several ships, syn- 
optic analyses and high resolution satellite imagery. The resultant 
examination yielded a description of an evaporation duct character- 
istic of the region as it was affected by frequently changing weather 
patterns. 

This study was on data collected from four research vessels 
eaepaeed im the JASIN experiment from mid-July to early September, 1978 
under open ocean conditions. The quantity of data and the large area 
of study permitted an examination to determine the nature of the duct 
and the atmospheric disturbances most likely to influence it. An 
important aspect of the experiment was the positioning of four ships in 
a triangle (180 x 180 x 200 km) during periods of intensive meteorologi- 
cal measurements. This constant spatial separation permitted accurate 
assessments of spatial as well as temporal variations over a large 
horizontal region as weather systems moved through the area. The data 
were used with expressions presented by Fairall et al [1978] to compute 


the associated oceanic evaporation duct height. 





II. BACKGROUND 


It has been long recognized that distributions of temperature and 
humidity in the lower atmosphere are responsible for anomalous electro- 
magnetic (EM) wave propagation [Katzin, 1947 and Kerr, 1951]  result- 
ing in either enhancement or degradation of Naval sensors and weapons 
systems. In particular, propagation between terminals reasonably close 
to the sea surface has been observed to exhibit significant variation 
of signal strength at ranges far in excess of the radar horizon. 

(Fig. 1) Such anomalous propagation effects were first seriously 
Studied during and immediately following World War II. In recent years, 
advances in mixed-layer theory have led to numerous computational 
approaches for real-time assessment of atmospheric refractive effects. 

Often relatively cool moist marine air extends vertically from 
the ocean surface to altitudes of a few hundred meters. Above these 
heights the air becomes much warmer and drier for a variety of reasons. 
Therefore, a transition laver exists in which the air temperature in- 
creases and moisture rapidly decreases with increasing height. The 
Pesultimemeradients can cause the EM refractive index to decrease 
rapidly enough to refract low incident angle radio waves downward 
relative to the earth and create ducting conditions. If the decrease 
in refractive index is great enough a surface-based duct occurs which 
will trap surface-originated radar signals and give greatly extended 


RoigesetOnedetecLion Of near Surface targets. 
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Figure la. Radar wave path under ''standard'' atmospheric conditions. 
Note path curves downward but at a rate less than the earth's curvature. 
Beyond-the-horizon target detection is not possible. 
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Figure 1b. Radar wave path under ducting conditions. Path curves down- 
ward at a rate exceeding the earth's curvature resulting in beyond-the- 
horizon target detection. 
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The air adjacent to the ocean is saturated and the relative 
humidity 1s near 100%, decreasing rapidly in the first few meters to 
ambient values which depend on varying meteorological conditions. The 
decrease of humidity causes the modified refractive index (M) to 
decrease with height initially. At greater heights, the humidity 
gradient decreases which will cause M to reach a minimum and thereafter 
increase with height (Fig. 2). The height at which M reaches a 
minimum value is called the evaporation duct height (Z*) and is a 


measure of the strength of the evaporation duct. 
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Figure 2. Relative humidity and modified refractivity M versus 
altitude for an evaporation duct. 
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The probability of occurrence of duct heights large enough to 
cause beyond-the-horizon detection capability for a particular radar 
varies according to geographic location, season and the time of day. 
The strongest ducts are generally found at tropical latitudes, during 


the summer season, and during daylight hours. 
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ITI. EVAPORATION DUCT MODEL 


Flux-profile relationships in the atmospheric surface layer 
(Z<50m) above the ocean surface have been ersten extensively and are 
well understood. Bulk methods can be used to determine surface flux 
scaling parameters from routine shipboard measurements. These scaling 
parameters can be applied to refractive equations to determine Z*. 
The method requires only mean meteorological conditions as input para- 
meters and is well suited to assessments using routine shipboard obser- 


Vaelons . 


A. EM PROPAGATION CONSIDERATIONS 

The refractivity N (in N units) of the atmosphere for EM waves is 
given by Battan [1959] 

N = 77.6 P/T + 3.73x10° e/T* co) 
where P is the atmospheric pressure, T is the absolute temperature 
and e is the partial pressure of water vapor. An expression for the 
vertical gradient in N is given by 

aids = Cl 6P/6Z + C2 éq/é6Z + C3 $1/d2Z rz) 
where q is the water vapor mixing ratio (q=.625e/P in gm/kg, P in mb, 
ieiniee, ang Cl=.5, €2=7.2 and C3=-1.35). A megative gradient in refrac- 
tivity at height Z causes a horizontally propagating EM wave to be bent 
toward the surface and if the gradient is large enough, the amount of 
bending exceeds the curvature of the earth and the wave becomes ducted. 
Since magnitudes of specific humidity gradients decrease with increas- 


ing height in the near-surface layer, the gradient is no longer strong 


f 





enough for ducting at some height Z*. This height is termed the evapo- 
rative duct height and is defined by the critical gradient necessary 
for EM trapping. 

dN/dZ = -.157 (m') 
using Equation (2) and noting dP/dZ= -oq (the hydrostatic balance for 
the atmosphere), Z* is then defined by the following equation. 

eto ewer? 6G/62 +°C5 61/62 (3) 
Since 6q/6Z and 6T/6Z are height dependent (approximately 1/Z) there 
exists a value of Z=Z* such that Equation (3) is satisfied. The bulk 
method allows one to relate 6q/6Z and 6T/5Z, which are needed to deter- 


mine Z*, to the surface layer fluxes and, hence, to mean parameters. 


B. BULK METHOD 

The theoretical framework for determining the evaporation duct 
height comes from knowledge of surface layer scaling parameters using 
the Monin-Obukhov similarity (MOS) theory. The theory relates the sur- 
face layer profiles of temperature, wind speed, water vapor and turbu- 
lence to the surface fluxes of momentum, sensible and latent heat. 
Vertical gradients of wind speed (U), potential temperature (T) or water 
vapor mixing ratio (q) are given by the following general expression 


6X/6Z = (X*/akZ) 9, (Z/L) (4) 


Where X = U, T or q, « = Von Karman's constant (0.35), a is a dif- 
fusivity constant, X* is the scaling parameter and 9(Z/L) is the MOS 
E@aptildtyecorrection. Ihe surface fluxes define the scaling parameters 


within the flux-profile relationships. Scaling parameters for 
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velocity, potential temperature and specific humidity are related to 
the normalized fluxes by: 
Ut = (-UTw')? is) 
(-T'w')/U* (°K) 
(-q'w')/U* (g/Kg) 


T* 


q* 
where u'w', T'w' and q'w' are the mean turbulent fluxes of momentum, heat 
and moisture respectively. A scaling length, the Monin-Obukhov length 
(L), 1s used to account for stability influences. Stability magnitudes 
are based on the height scaled by L, i.e. Z/L. Negative and positive 
Z/L values corresponded to unstable and stable conditions. L is defined 
by the scaling parameters as 

L = TU* /(<g(T* + .61q*T)) (5) 
The bulk method is based on integral forms of Equation (4). Integrating 


Equation (4) from Z = to some reference height, Z, an expression 


on 
for the value of a parameter at the measurement height is obtained 
(Businger, 1973]. 

X(2) = X(Zy) + X*/ax(In 2/2, - (2/1) (6) 
where (Z/L) is the profile stability function and subscript 

X = U, T, or q. (Note: a and ~ depend on X) 

Utilizing Equation (6), a method of estimating surfaces fluxes in- 
volves the scaling parameters and roughness lengths, 2x" Appropriate 
values of temperature, relative humidity and wind speed are obtained 


at both the sea surface and at some reference height, Z. Equation (6) 


can be used to determine scaling as a function of the differences. 


it 





U* 


(cU) / (In (2/2) - vi (2/1) 
(T-Ts) ak/ (In(Z/Z.5) -be (Z/L)) 
(q-qs) ak/ (1n(Z/Z5n) “v2 (2/L)) 


a was set equal to 1.55 for T and q from Businger et al [1971]. The 


T* 


ol 


quantities ox and Zor are roughness lengths for velocity and tempera- 
ture profiles. It is convenient to rewrite these equations in the 


Beweralized drag coefficient form 


1 
“2 


X* = Cy? (X(Z) - X(0)) @ 
The neutral stability drag coefficient can be defined as 
L. 
Coy = ak/In(Z/Z,y) (8) 
liens the drag coefficient becomes 
1. 
Meier (i ((C..") fox) W (2/L)))2 (9) 
ite thas form, 2 Ox is related to the neutral drag coefficient by 
eh 
Cy = ak/  In(Z/Z,.) (10) 
where roughness length 
2 
ae x ON) (11) 


Equation (5) can be rewritten to give the MOS stability parameter, 





wi. as 
tee 4) / Z/L)) )2 
ee ee (12) 
U/L = (2/1), : 
(1 = (Cp) fax (eo (Z/L))) 
where 
Cc. T(Z)-T(0)) + .18(q(Z)-a(0))) 
cat CIN E(T)-TEO)) + -18((Z)-a -_ 
(2/L) = a U(Z) 
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Note that (2/1) , 1S an initial estimate of stability in terms of 
differences and neutral drag coefficient. Final values of Z/L and 
scaling parameters can be found from a simple iteration process. 
Stability correction functions used were: 


Wellocley profile: 


eee) eee ime + X)/2) + In ((1 4 x )/2) 
-2 ARCTAN (X) * 11/2 Z2/L<0 

X= (1 -15 Z/L)* 

Wi (Z/L) = -4.7(Z/L) Z/L>0 

Temperature profile: 

Uae = 2 ine ((l + X) /2) Z/L<0 
X= (1 -92/L)%. 

W2(Z/L) = -6.5(2/L) Z/L>0 


The velocity drag coefficient, CN in a near-neutral atmosphere 
was Calculated from the wind-speed dependent formulation given by 
Kondo [1975 ]. A typical value of Cy at Z=10 m was 1.3x10 ? m/s 
which yielded 2 =6x10 * m. Wind speeds and drag coefficients are 


given in Table [. 


Ay 





TABLE I 


CoN VERSUS WIND SPEED AT 10m[ KONDO, 1975 ] 





3 
U(m/s) Coy x 10 
Mood ac ieOSxUe. 2 

5:20 jet eavioos x U 

5.0-8.0 Sieur 2c U 

8.0-25.0 lee elo x 
The temperature and moisture drag coefficients, CoN and CON? ee = POM 
were estimated as 1.3x10 ?[Davidson et al, 1978]. For Gp=1. 35 


[Businger, et al, 1971], a roughness length of Z_.=2.0x10 °* m, was 


Ol 
obtained, which is assumed to be independent of wind speed. 


C. Z* EQUATION IN MOS FORM 

In the surface layer, the gradient of N can be written in terms of 
and q* as 

ahydZ = -.032 = 9(Z/L)/(akzZ) (7.24% - 1.3T*) (14) 


which arises from Equation (3). The duct height, 2*, is obtained from 


Equation (14) by setting dN/dZ = -.157, the critical gradient for ducting. 
igecg: = 1.50%) 
7% = a a 7 Ty) (15) 
@ic( 125} 


Given q*, T* and L, Equation (15) was solved iteratively to obtain 2*. 
No ducting occurs unless 


7.2q* -1.3T*< 0 
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IV. DATA ACQUISITION 


The data used in this analysis were unique in that they represented 
accurate measurements of the atmospheric surface boundary layer on a 
horizontal scale of 1-200 km. All data were gathered from scientific 
research vessels whose primary cruise objectives were to observe, dis- 
tinguish and quantify fluxes and mixing in atmospheric and oceanic 
boundary layers. The following are descriptions of the measurement 
platforms, their locations, the data accuracy and the intercomparison 


of data between ships. 


A. PLATFORMS AND LOCATION 


Data for this analysis were taken during the JASIN experiment in the 


— _— —_—- 


summer of 1978. Data from four eiits, invamed ae experiment were 
selected on the basis of two objectives: 

1. Examine as much of the period as possible. 

2. Examine variation over nearly fixed spatial separation. 
The four ships selected were the METEOR, the JOHN MURRAY, the GARDLINE 
ENDURER and the HECLA. Figure 3 illustrates the profile and sensor loca- 
tions for each ship. Hourly observations were made of air temperature 
(T), wet bulb temperature (Tw), sea (bucket) surface temperature (Ts) 
and wind speed (U). The triangle configuration of the ship locations, 


depicted in Figure 4, provided quite good areal coverage. The ships 


were positioned approximately 200 km apart at each corner of the 


i 
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Figure 3. Ship profiles and sensor locations: (A)- JU, 
(B= Ti lw ande(e je is. 
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Figure 4. The JASIN area in relation to the British Isles. 
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triangle and one ship was located at the center. Nominal positions 
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B. DATA COLLECTION 

The subject observations were recorded on log sheets and plots of 
three-hourly values were prepared by JASIN participants. The plots 
from the GARDLINE ENDURER, JOHN MURRAY and HECLA, were the primary 
source of data used in this study. Air temperature, wet bulb tempera- 
ture, sea surface temperature were accurate to 0.2 mee PRessure= CO 
1 mb and wind speed to 1 kt. Additionally, data used from the METEOR 


were available in computer printout. 


C. DATA QUALITY AND INTERCOMPARISON 

During the period 15-22 July, on 19 and 20 August, and on 
5 September, the four ships were positioned for an intercomparison of 
meteorological measurements. Hourly (half-hourly during formal inter- 
comparisons) values of wind speed, atmospheric pressure, dry and wet 
bulb air temperatures and sea surface temperature were compared between 
vessels approximately one kilometer apart. These data were initially 
evaluated for coding errors and inconsistencies resulting from position 
changes. Mean and standard deviations of each variable were calculated 
for each participating vessel. Time-series of inter-platform sensor 


disagreements were produced from mathematical differences between 


oe 


a 
= 





observed variables. Interplatform means and standard deviations of 
these differences were calculated and, assuming the disagreements to be 
normally distributed, the 90% confidence interval of the estimates of 
differences were calculated. 

Results derived from this procedure were time records of inter- 
platform sensor differences. Constant disagreement throughout JASIN 
implied the sensor performed in a consistent manner and it remained 
only to explain these errors in terms of known disparities, such as 
non-uniform sensor heights or instrument bias. A variation of the 
interplatform corrections required further analyses either in terms of 
atmospheric stability or in terms of trends of sensor drifts. The 
latter were established by regression techniques. Inconsistencies 
unexplained or uncorrected by these considerations remained as errors 
in the four-ship system. 

The METEOR Fata were selected as nner to Sich date Senn the 
other vessels were adjusted. However, the METEOR's sea surface tempera- 
ture values were corrected with data from the other three ships. The 
METEOR was chosen because it was involved in more intercomparisons, 
because it had a preferred method for station-keeping, and because of 
its proximity to other JASIN ships and buoys. The corrections accord- 
ing to measured parameters and individual ships are summarized below 
{Macklin and Guymer, 1980]. 

1. Wind speed (m/s) 

METEOR = METEOR 

HECLA = 0.89 (HECLA) + 0.3 

MURRAY = MURRAY 

ENDURER = ENDURER 

OR 





2. Pressure (mb) 


METEOR = METEOR 
HECLA = 


HECLA + .40 From 0900, 


HECLA + .08 To 0600, July 28 


July 28 


MURRAY = 1.00640 (MURRAY) - 6.23 


ENDURER 


3. Dry Bulb Temperature (°C) 
METEOR = METEOR 


HECLA = HECLA - 0.3 
MURRAY = MURRAY - 0.5 
ENDURER = ENDURER - 0.5 


4, Wet Bulb Temperature (°C) 


METEOR = METEOR 

megan = HECLA - 0.3 
MURRAY = MURRAY - 0.3 
ENDURER = ENDURER - 0.1 


5. Sea Bucket Temperature (°C) 


ENDURER + .38 To 1700, August 19 
ENDURER + .38 After 1700, August 19 (P less than 1000) 
ENDURER + .48 After 1700, August 19 (Pp greater than 1000) 


PeteOR = METEOR + 0.3 To 1300, July 27 
= METEOR FOC; sJiive27 to Ol00Fenueust 3 
= METEOR - 0.4 0100, August 3 to 1200, August 9 
= METEOR From 1200, August 9 

HECLA = HECLA 


MURRAY = MURRAY 
ENDURER = ENDURER 
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V. ANALYSIS PROCEDURES 


Data analysis procedures were designed to determine 1) the repre- 
sentative Z* values for the summer regime, 2) meteorological factors 
which most affected Z*, 3) the frequency of occurrence and duration of 
ducts, and 4) the extent of horizontal homogeneity of the duct. To 
achieve this, the data were first corrected, then used as input to 
bulk formulae to compute Z*. 

A data file was constructed of the most representative Z* values 
and corrected Ts-T, U and RH at three-hour intervals. Winds were 
reduced from anemometer height to temperature sensor heights using the 
Stability-corrected logarithmic wind profile [Businger et al, 1971]. 
The periods examined and sensor heights are summarized in Table II. 
(Dates are in year-month-day-hour format) 


TABLE II 
SHIP DATA PERIODS AND MEASUREMENT HEIGHTS 





Sill P STARTING ENDING ANEMOMETER HGT (m) Z({m) 
MURRAY foUPa SUZ 78072103 11 6 
METEOR 78072106 78080912 oe 11 
MURRAY 78080915 78081012 sla 6 
METEOR 78081709 78090506 ae i 
ENDURER 78090509 78090606 15 + 


Statistical and time-series analyses of Z* were performed along 


with subjective analyses of influencing meteorological conditions. The 


latter were performed with individual observations, synoptic weather 


maps, and satellite imagery (visual and infrared). 


were correlated to the time-series of Z* values. 


a 


These descriptions 


Joint probability 





density distribution and conditional mean function results were chosen 
to examine correlations of at least two weather parameters to Z* 
values. Finally, comparisons were made of Z* and all atmospheric 
measurements to assess the horizontal homogeneity of the duct for 


periods when ships were spatially separated and at fixed locations. 


VI. PRESENTATION AND INTERPRETATION OF RESULTS 


A. GENERAL ANALYSIS OF SUMMER PERIOD (13 July - 6 September) 

Figure 5 1S a time-series comparison of evaporation duct height 
(Z*), air-sea temperature difference (Ts-T), wind speed (U), and 
relative humidity (RH) for the period of study, 13 July to 6 September, 
with the exception of the interval between 1200 GMT, 10 August and 0900 
GMT, 17 August. This time series was formulated from at least one ship 
located in the JASIN array and provided a near-continuous analysis of 
the evaporation duct height (Z*) at three-hour intervals. General 
conditions during the period will be described in the following sections. 

1. Synoptic Conditions in JASIN Area 

Daily weather maps for the Northeast Atlantic are shown in 
Appendix A. Julian day and date (day of the month) are shown on each 
map (see also Table A-I). Studies of the daily surface pressure charts 
and hourly observations [ Guymer and Taylor, 1978] showed that there 
were essentially four synoptic situations affecting the JASIN area 


during the entire period. 
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a. Moist Anticyclonic (A) 

The periods of 14-18 July and 22-30 August were affected by 
moist westerlies or northwesterlies on the northem flank of an anti- 
cyclone to the west of Ireland. Weak frontal waves were embedded in 
this flow and occasional incursions of drier air took place. Winds 
varied between 3 and 12 m/s. 

b. Mobile Westerly Conditions (B) 

The periods of 19-21 July and 30 August - 2 September were 
characterized by depressions moving through the area bringing cold and 
dry air. The disturbances were very weak. 

c. Cyclonic Southwesterly (C) 

During the periods 22-29 July; 11-21 August and 2-5 Sep- 
tember the depressions became slow moving to the south of Iceland and 
muomes or troughs from the southwest affected the region. Generally, 
the fronts were occluded and slow moving by the time they reached the 
area, with only moderate winds. However, winds which were the strong- 
est experienced during the entire JASIN period reached speeds in excess 
meee son 22/26 July, 18/20 August and 5 September. The air was 
Significantly warmer than the sea. 

d. Blocked (D) 

During the period 30 July - 10 August a blocking low formed 
over northwest France and moved very slowly northeast; to the west a 
ridge persisted. Low level northeasterlies were colder than the sea 
surface but convection was inhibited by warm air aloft associated with 
weak, slow moving occlusions which also produced occasional drizzle. 
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Towards the end of the period, the ridge intensified and moved eastward. 
e. Special Events 
some interesting atmospheric occurrences were noted as 
events (refer to events marked on Figure 5) and are listed below. 


1. Relatively cold, dry air (14, 20 July; overnight 1-2 August; 
fone co, oi August). 


2. Significant moisture flux in near-neutral conditions (30 July; 
23-25 August; 1 September). 


=> 


5. A warm frontal wave on 25 August and a cold front preceded by 
convective instability on 28 August. Surface fluxes were very small. 


fee aids in| excess of 15 m/s (22 July; 18, 19 August). 
ei ewmsions Of Warm, dry air (14, 18, 25, 28, 29 July). 
The following are descriptions of events which caused warm, dry air 
(not a normal occurrence over the open ocean) to occur. 
(1) Afternoon Heating. On 14 July, both the MURRAY and 
HECLA reported increasing:air temperature but MURRAY reported three 
consecutive wet bulb temperature readings markedly lower than HECLA. 
On 18 July, MURRAY and METEOR had similar air temperature readings 
while HECLA was 1.4°C lower. All three measured similar values for 
wet bulb temperature. On both occasions, winds were light under 
prevailing high pressure. 
(2) Subsidence and Continental Trajectory. On 25 July, 
a weak ridge moved over the area as a low west of Ireland moved north- 
east. Winds backed from the northwest to the east-southeast. Warm,dry 
air, a result of subsidence from earlier ridging, flowed from the 


British Isles and influenced the area for six hours before the approaching 
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occlusion became dominant. 

(5) Short Wave Troughs. Throughout 28 and 29 July, a long 
wave trough was oriented north-south from Iceland. Weak short wave 
impulses on both days provided surges of dry air, each lasting nine hours, 


while generally maintaining warm southerly flow. 


2. Sea Surface Temperature (TS) 

surface water characteristics were associated with the North 
Atlantic Drift, an extension of the warm waters of the Gulf Stream. The 
mean water temperature in the vicinity of the JASIN triangle was 12.2°C. 
The general areal distribution of the surface temperature for the 
sampled periods is depicted in Figure 6 [Liu §& Katsaros, 1978]. 
The Northeast corner of the triangle was always the coldest position 
(11.2°C). The Northwest and South corners of the triangle were always 
the warmest (between 12 and 13°C). A large tongue of cool water 
(11.8°C) from the north, intruded south and south-southeast through 
the center of the triangle. It persisted for the entire period result- 
ing in minimal temperature variations at the center position (about 
0.2°C). The most pronounced Ts gradients occurred along a line between 


the NW and NE stations. 


3. Representative Z* Values 
A statistical analysis of observed parameters and calculated 
values from 714 independent observations was performed to determine a 


representative Z* value (Table III). The mean evaporation duct height 
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was 4.2 meters with a standard deviation of 3.9 meters. This height 
appears low when compared with the climatology of the region. Sweet 
[1979] determined an’ evaporation duct climatology at Ocean Station 
ALFA (62N/33W), 1100 km upstream of the JASIN array, and Ocean 
Station INDIA (58N/19W), 400 km southwest of the JASIN array. The 
following statistics are a comparison of the JASIN evaporation duct, 
using the method previously described and climatology which used bulk 
aerodynamic formulae and a method originated by Jeske [1971] and modi- 


fied by Hitney [1975 ]. 





ALFA JASIN INDIA 
MEAN (m) 5.8 ee 7 
MEDIAN (m) 5.7 3.8 7 

TOR (m) 3.4-8.3 Oe Ser 4.1-9.8 
# OBS 7058 714 6758 


The interquartile range (IQR) represents the lower and upper bounds for 
the middle duct heights (25% of the ducts have heights below the lower 
number and 75% have heights below the upper number) and can be viewed 
as an indicator of the spread or variation of the duct height. 
TABLE III 
Sie mioiiCs FOR JASIN PERIOD 1s JULY - 6 SEP 1978 
(714 OBSERVATIONS) 


is T TS-T U RH Z/L T*® U*  q* ~ 28 
ese) er Oe (m/s) ) (°C) (m/s) (g/kg) () 


MEAN ec lise .5 Oona U eee) = UZ 


SIGMA (Owes: | A ae 8 0 ee 15 sOSMFZeo 
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4, Atmospheric Factors Affecting Z* 

Since Z* is a computational artifact determined by four meteor- 
ological parameters, Z* (Ts-T, Rh, U, Z/L), multivariable analyses 
were used to determine the degree to which individual variables affected 
duct strength. Joint probability density distributions were computed 
for each of the parametric interrelationships to examine the pattern of 
atmospheric factors as they are related to the occurrence of Z* values. 

Based on probability theory [Batchelor, 1953] and utilizing a 
technique described by Holland [1968] , the distribution of Z* was ana- 
lyzed statistically against the distribution occurrence of a combination 
of meteorological variable such as Ts-T versus U, (Figures 7-12). 
Individual moments and distributions for each variable were also computed 
in the analysis. The following data were tabulated: 
Mean 
Variance 
Standard Deviation (Sigma) 
Skewness 
Kurtosis 
Dimensional Class Limits by Sigma/4 Intervals From -4.0 Sigma 
to +4.0 Sigma 


7. Density of Observations for Each Occurrence by Sigma/4 Class 
iieerva | 


HNN GN Fe 


Individual class and joint class intervals were computed for 
each observation set. Two indices were computed for each data pair 
representing its position in an array (18 x 18) of joint class inter- 
vals (Sigma/2 x Sigma/2) centered on the mean values of the two 
variables. Class intervals were determined by dividing the standard 
value of one of the joint variables by 0.5. Pairs of independent 


variables were then used to compute joint probability density functions 
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Figure 7. Conditional mean (plotted numbers) of Z* for 
joint-frequency distribution (delineated with solid labelled 
isolines) of wind speed, (U) and air-sea temperature differ- 
ence, (Ts-T). Integers along the right and bottom edges are 
sums of occurrences in respective rows or columns. The total 
sum is 714 but may be less if the compared variables are 
truncated when out of range. 
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Figure 8 (top). JDF of Z/L and Figure 9 (bottom) JDF 
of RH and Ts-T to CMF of Z*. 
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Figure 10. JDF of wind speed and relative humidity 
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Figure 11. JDF of relative humidity and stability 
SOUCMEROTa2s: 


36 





e 1302 206 e ° ° e e e e e ° 






e 40.0 —— tel = ° e 
6.2 tea? aac 305. RO iso 13) 40 BP 40> 5-0 4.0 p ayaa 609 e 
308 308 5.0 ae 9 9 e e 

10 


109 3e3 508 e e °o 3S o e e e e e 





105 ® e 9 e e e e e @ e Go ® 
e 203 e e Qo 9 e e e e e S ® e 
e 205 e e o e e e e e @ e e e 
-3.1 
e 2e3 tC e e e e e e e e @ e e 


18 49 107 149 491 96 76 50 19 7 5 2 


0.8 37 6.6 95 124 15.3 187 


U (m/s ) 


Figure 12. JDF wind speed and stability to CMF of 2Z*. 
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(JDF). Similarly, a conditional mean function (CMF), defined as the 
average values of Z* for all observation times when the two independent 
variables had values within the joint intervals, was defined for each 
joint class interval. A subjective analysis of the JDF's and the CMF's 
was performed to distinguish relationships as being from computational 
artifacts or from accidental occurrences. 

pear -oea lemperatume Difference {Ts-1T): (Figures 7, 8 and 9) 

The air temperature averaged 0.5°C colder than the sea sur- 

face temperature. Z*(Ts) calculations were largely influenced by the 
observing ship's location in the JASIN area (Fig. 6). Air temperature 
variations were primarily influenced by migratory weather systems. The 
most notable positive deviation of Z* occurred when Ts-T was less than 
zero (-.5°C), with low relative humidity (RH<85%) and stable atmospher- 
ic conditions (Z/[>0). However, under these conditions duct strength 
was extremely sensitive to humidity and stability changes. If RH in- 
creased Sigma(RH) or Z/L decreased Sigma (Z/L) then the evaporation duct 
virtually vanished. Generally, Ts-T had little effect on Z* when 
RH>90%. As Ts-T approached zero or became positive, more positive 2Z* 
deviations occurred for a wider range of RH e.g. (50-90%). The poorest 
ducting conditions occurred when T exceeded Ts more than 0.5°C and RH 
was in the high 90%'s. In near-neutral conditions, as Ts-T increased 
Z* tended to increase fractionally. Generally then, Ts-T was more of a 
determinant than stability, except in the special situation of warn, 


moist alr over cooler water. 
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b. Relative Humidity: (Figures 9, 10 and 11) 

Average humidity was 90% and ranged from 49% to 100%. 
Significant changes in humidity were directly due to the passage of 
large-scale weather systems. As relative humidity dropped to 803, 
marked increases in Z* occurred. The largest Z* values occurred when 
dry (RH<80%¢) and warm air overlay cool water (Ts-T <-.5°C). For the 
relative humidity range (50-90%), larger than average Z* values 
occurred when the air was 1.0°C cooler than the sea surface. Lowest 
ducts occurred in warm, moist meteorological conditions (RH:85-100%, 
ie=1<-Z2 to 0). 

c. Wind Speed (U): (Figures 7, 10 and 12) 

Wind was most influential on Z* values in low wind situa- 
tions. In light winds, changes in atmospheric stability (Z/L) drasti- 
cally changed the duct. As seen in Figure 12, a change in Z/L by 
Sigma (Z/L) changed Z* from 3 to 40 meters for winds near 5 m/s. In 
moist air (RH>90%), Z* was low regardless of wind speed. In dryer air 
(RH<85%3) wind caused increased Z* as U increased from 5 to 12 m/s. 

demmocapllity (2/L); (Figures 8, Jl and 12) 

Average stability was near-neutral. The greatest deviations 
of Z* occurred when Ts-T<0 (warm air) and light wind (U<5m/s). Other- 
wise changes in Z/L, hence changes in Z*, were small as winds increased 
or colder air was advected over relatively warmer water (Ts-T>0). 

Figure 13 illustrates the interrelationships that the 


combination of selected parameters had on Z*. The ranges of Z* are 


a9 








depicted for the primary and secondary extremes encountered in the 
joint probability density distribution. With the exception of the 
dependent parametric interrelationship of stability (Z/L) and air-sea 
temperature (Ts-T), all parameters caused large values of Z*. The 
primary extreme range was a reflection of the special situation in 
which large values of Z* were generated due to atmospheric conditions 
where winds were light and warm, dry, stable air overlay relatively 
cooler waters. A more significant indication of atmospheric influence 
was the secondary range. This range reflected more accurately the 
effects due to characteristic changes in atmospheric variables with 
synoptic conditions of a summer regime. 

The most influential meteorological factors were air-sea 
temperature and relative humidity, followed closely by the interrela- 
tionship of wind speed to air-sea temperature. Because the secondary 
range of wind speed to relative humidity was small, this indicated 
that large positive Z* deviations were influenced most significantly 
by air-sea temperature. Stability only became a significant factor 
in light winds. Table IV summarizes the relationship of duct strength 


to atmospheric parameters. 
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Figure 13. Parametric Interrelationship to Duct Height 
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TABLE IV 


DUCT STRENGTH VS. ATMOSPHERIC PARAMETERS 


z% RH (3) Dene Tsar C) Z/L 
STRONGEST < 80 < 5 nS > 0 
STRONG < 85 Baiky > 0 
WEAK 85-100 -2 to 0 > 0 
WEAK > 90 > 0 


5. Z* OCCURRENCE AND DURATION 


—_——_——- 


peter rans of duct occurrence and duration of duct heights are 
shown im Figure 14. ‘This figure indicates the occurrence of strong 
evaporation ducts during the summer is small and short-lived. Duct 
heights of 6m or less occurred 78% of the time. Rarely did significant 
ducts occur for more than several hours. Figure 5 shows the changing 
weather conditions throughout the summer and the effect on duct height. 
Rarely did one weather pattern dominate or persist. Consequently, the 
occurrence of strong ducting was irregular. It was most closely asso- 
Ciated with atmospheric conditions that caused either cold, dry air or 
warm, dry air to move over the area. Incursions of relatively cold, 
dry air, usually from the northwest (Figures 21-25, 28-33 and 40-45), 
mised rapid formation of evaporation ducts (9-12m) which rarely lasted 
longer than 9-12 hours. Dry air did not last long over the area before 
it was modified. Incursions of warm dry air were evident and produced 
brief (3-6 hrs) but strong evaporation ducts (>8m). Conditions less 


favorable for the formation of ducts were associated with cyclonic 


southwesterly flow of warm, moist air (Figure 20). 
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Figure 14. Frequency of duct occurence and duration of duct 
heights. 
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Figure 15 and 16 compare the distribution frequency of Z*. Of the 
three locations, JASIN indicates a higher percentage of low duct 
heights, Ocean Station INDIA shows its mode (group that occurs most 
frequently) in the 6-9 m range, the highest of the three. Additionally, 
INDIA had a greater percentage of Z* in the 9-18m range. Ocean Station 
ALFA located off the east coast of Greenland and in an oceanic region 
Pe confluence where the cold Greenland current meets the cool Irminger 
current, indicated a mode of 3-6m and a generally higher percentage of 
Z* greater than JASIN but less than INDIA. 

This difference between the evaporation duct climatology in the 
three areas appears to have been due to (1) the proximity of ALFA to 
the Greenland landmass and the associated, extremely cold water and 
(2) INDIA's location in the deeper and warmer water of the North 
Atlantic current. Interpretation of visual and infrared satellite 
imagery indicates (Figures 28-33) the flow of cold, dry air off 
Greenland, coupled with cold sea surface temperatures, causes higher 


iE 


values at ALFA. due to the large difference in surface humidity 

and atmospheric humidity. As the dry air mass moves across the North 
Atlantic atmospheric humidity changes are relatively small as it 
approaches both the JASIN and INDIA locations. The key difference 
between the two locations is Ts. The warmer water at INDIA causes 
— surface humidity, thus a larger difference in atmospheric and 
surface humidity resulting in the occurrence of stronger ducts. JASIN, 


because the air mass has become more moist since leaving ALFA, but 


enters an area of relatively warmer water, experiences the smallest 
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Figure 16. Frequency of 2* occurrence (climatology) at 
North Atlantic Ocean Stations nearest JASIN. 
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Ocean Stations in the North Atlantic: ALFA 


(62N, 33W) and INDIA (58N, 19W). 
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difference in atmospheric and surface humidities, therefore the weakest 


@uets of the three locations. 


B. HORIZONTAL NATURE OF THE EVAPORATION DUCT 

The three periods (7-8, 23-25 August and 30 August - 1 September) 
examined corresponded to times when ships were positioned at all corners 
of the JASIN triangle. The objectives of analysis for these periods 
mene CO: 

1. Determine the horizontal homogeneity of the evaporation duct. 

2. Determine which meteorological factors enhanced or degraded areal 

ducting. 
The analysis procedure for this was to: 

1. Compare computed Z* values at different ship positions and con- 

comitant meteorological parameters to assess similarity. 

Z. Analyze synoptic meteorological conditions utilizing satellite 
imagery (visual/infrared) , ship observations and synoptic weather 
maps. 

3. Analyze the temporal variations or consistency of horizontally 
homogeneous conditions. 

4. Relate the above fields to the spatial nature of the evaporation 

auet. 

Area homogeneity was determined by comparing Z*, Ts-T, RH and U 
for ships positioned at each corner of the triangle. Azimuthal sectors 
were defined as the legs of the triangle (NW sector: NW-S, NE sector: 
NE-S, E sector: NE-NW). Using the South (S) position as the reference, 
differences were determined relative to Northwest (NW) and Northeast 
(NE) locations (NW and NE sectors, respectively). Similarly, using the 
Northwest observer as reference, differences were determined relative 


to the Northeast location (East sector). Differences between locations 


(ships) relative to sectors are tabulated in Tables V through VII. 
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1. 2* Range in Homogeneous Atmosphere 


A necessary step to interpret the results of the tabulated 
differences is the performance of an error analysis to determine the 
possible range of computed Z* which could occur in a homogeneous 
atmospheric situation due to measurement inaccuracies. An assessment 
was required to ascertain the maximum deviation allowable in Z* for the 
three-ship system under nearly homogeneous atmospheric conditions. 
Additionally, a sensitivity analysis of Z* values was performed for 
each inherent measurement uncertainty. 

Previous studies on bulk measurements [Davidson et al, 1978] 
have concluded that representative inaccuracies in bulk parameter 
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Five observation times were chosen to perform the error analysis. 
These corresponded to times when air temperature, relative humidity and 
wind speed were nearly similar for all ships stationed on the triangle. 
Benchmark Z* values were computed from each of these representative 
homogeneous atmospheric situations. Atmospheric conditions at 1500 GMT 
August 23, 24 and 0600 GMI 25 August were similar at all three ship 
locations. Mean values of bulk parameters were used to compute a Z* 
at each corner of the triangle assuming no measurement inaccuracies. 
A similar analysis was performed on data taken at 0600 and 1500 GMT 31 


August. The following gives the results of the two periods: 
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CASE I 





is = 12 °C eee. TC RH = 90% U = 8.1 m/s 
SHIP Z(m ae Ee Gue 

ENDURER 4 =. 0055 A a4 

MURRAY 6 Se oe Ao 

METEOR el 20159 Aes 

CASE II 

is =a e5r"C T = isc RH = 75% U = 10 m/s 
SHIP Z(m) eis 2* (m) 

ENDURER 4 Sen ees Gea 

MURRAY 6 02.5 jib) call 

METEOR oe -,045 9.7 


The Z* values computed in each of the above cases represented 
reference Z* values for each measurement height, Z (e.g. 4, 6, 11m) in 
a spatially homogeneous, mixed layer. Measurement uncertainties listed 
below were introduced into the Z* calculations and the maximum error 
computed. Percentage values for both cases reflect the change in Z* 


due only to the measurement errors of each individual parameter. 


CASE I CASE II _ MEASUREMENT ERROR (ms) 
U 7% 33 10% 
RH 25% 89 3% 
Ts 14% 2% 6 


The cumulative error in the worst case situation, where measurement 
uncertainties were between two ships in a horizontally homogeneous 


atmosphere, produced the following difference in 2* values: 
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Therefore, it was concluded that if the atmosphere were homo- 
geneous and the only uncertainties were in measurements of T, Ts, RH 
and U, then the greatest difference in Z* that any two ships would 
have experienced was 2.9 m. 

2. Case Studies of Horizontal Homogeneity 

Three periods were examined, in view of the error analysis, when 
ships were stationed at the corners of the JASIN triangle (approximately 
200 kn apart). Simultaneous meteorological observations, taken at three- 
hour intervals, and their respective Z* values were compared. In each 
case, an average or better than average evaporation duct existed. Exten- 
Sive periods of below average ducting were not considered. 

The examination included consideration of the general synoptic 
conditions, differences in local meteorological conditions between ships 
(and the factors which most influenced Z*), the extent of the duct over 
the entire triangle (areal homogeneity) and between corners (sector homo- 
geneity). The maximum Z* difference (2.9m) between ships (due to 
measurement errors as previously described) was the criterion for duct 
homogeneity. Therefore, if the maximum difference in 2* between all 
three ships was less than or equal to 2.9m, areal homogeneity was 
assumed to exist. Similarly, if the criterion was valid between corners 


of the triangie, sector homogeneity was assumed to exist. 
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a. Case [I 
(1) Period. 7 Aug (0000 GMT) - 8 Aug (2100 GMT) 
(2) Observer Positions on Triangle. 
Northwest (NW) - ENDURER 


Northeast (NE) - MURRAY 
south (S) - METEOR 


(3) Synoptic Conditions - Blocked. A ridge persisted over 
the area as it intensified and moved Slowly east (Figure 19). Northerly 
surface winds were moderate in the eastern region and lighter in the 
west. Low level instability prevailed as cold air was advected over 
warmer water. However, convection was inhibited by warm air aloft as 
an occlusion moved slowly over the area (Figures 25-27). 

(4) Results of Tabulated Differences - (Table V). Obser- 
vations at NW and NE were more humid than S. Winds were stronger at 
NE than NW and S. Air-sea temperature differences were greatest at S. 
The duct heights on the average were 2m higher in the south than the 
two northern positions (i.e. 6.9m vs. 4.9m). 

(S) Areal Homogeneity. Similar evaporation duct strengths 
occurred simultaneously at the three stations 30% of the time with the 
most sustained conditions lasting 12 hours. Figure 18 illustrates by 
sectors the percentage of homogeneous ducts and the longest duration 
experienced. Northerly flow coincided with similarity between the two 
northern positions while the greatest disparity occurred between 
southern and western positions. The most persistent horizontal homo- 
geneity occurred between NW and NE where humidities and air-sea 


temperature differences were most similar. Poorest conditions existed 
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between NW and S due primarily to the dissimilarity in Ts-T and RH and 
the consequent difference in air-sea humidity (Figures 44-47). 
De Case II 
(1) Period: 22 August (1500 GMT) - 25 August (1500 GMT) 


(2) Observer Positions on Triangle: 
Northwest (NW) - ENDURER 


Northeast (NE) - HECLA 
south (S) = MEMEO 

(3) Synoptic Conditions - Moist Anticyclonic. An advanc- 
ing cold front passed over the area the afternoon of 23 August (Fig- 
ures 21, 28, 31). Cold, dry northwesterly air immediately followed the 
frontal passage and prevailed through 24 August (Figures 22,32, 33). 
Humidity gradually increased from the west over the area by 25 August 
as a warm front approached (Figures 23, 34, 35). 

(4) Results of Tabulated Differences -(Table VI}. Generally, 
the greatest differences experienced were between NW and NE stations 
due to the west to east movement of the air mass over the area. Winds 
and humidity were similar at all observing stations. Air-sea differ- 
ences were least at the NE position where water was the coldest, hence 
weakest ducts were generated. Duct heights ranged 1.1-12.2m and 
averaged 5m over the area. 

(5) Areal Homogeneity. The area experienced homogeneous 
ducting 44% of the time with the longest duration 9 hours. Homogeneity 
was nearly equal for all sectors with the northeast experiencing 


slightly stronger ducting and lasting for 24 hours following frontal 


passage (Figure 18). The greatest 2* difference across the front was 
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5.8m i.e. 6.6m (NW) vs. 2.8m (NE) Horizontally homogeneous ducting 
lasted for nearly 15 hours following frontal passage before another 
impulse of dry air entered the western area. It took 15 hours for the 
dry air to move across the JASIN area and strengthen the duct in the 
east (Figures 48-55). 
Sopecase i]t 
(1) Period. 30 August (0000 GMT) - 1 September (2100 GMT) 


(2) Observer Positions on Triangle. 





NORTHWEST (NW) - ENDURER 
NORTHEAST (NE) - HECLA 
SOUTH (S) - METEOR 
(3) Synoptic Conditions - Mobile Westerly. A weak 
depression moved over the region (Figure 19). The JASIN area was in a 
warm sector on 30 August as a cold front passed early on 31 August. A 
warm front approached from west larenen 1 September (Figures 40-43.2). 
(4) Results of Tabulated Difference - (Table VII). The 
average duct height was 6m. The most intense ducts occurred on 31 
August after the cold front passed. They occurred in the west and 
south where sea surface temperatures were highest and air-sea tempera- 
ture- differences were greatest. On the average, relative humidity 
was lower to the northeast than the western and southern region which 
were largely under the influence of warm sector of the low. 
(5) Areal Homogeneity. Due to the high mobility of the 
weather systems, areal horizontal homogeneous ducting occurred only 29% 
of che time with a maximum duration of 9 hours during very low ducts. 


Sector homogeneity was poorest between the northeast and northwest 
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positions and best between south and northeast (Figure 18). The maxi- 
mum difference in Z* across the cold front was 9.5m i.e. 11.7m (NW) 
vs. 2.2m(S). The strongest duct heights associated with any weather 
system were generated on 31 August following the cold front. Though 
outside the limits of horizontal homogeneity criterion, 2.9m, duct 
heights of 8.6-13.4m occurred at all stations for 12 hours before 
weakening in the east (Figures 54-59). 

5. Summary of Horizontal Variation of Z* 

The periods of the three cases studied were chosen because 
higher than normal ducts occurred for reasonably long periods. The 
horizontal nature of significant ducting in the region could be assessed. 
Each case was unique in that strong ducting was the result of differing 
Synoptic meteorological conditions. 

The strongest ducts occurred as mobile westerly depressions 
moved through the area; however, horizontal homogeneity was poorest. 
Air-sea temperature differences and relative humidity were the most 
important parameters affecting duct strength. Larger Z* values were 
experienced at the NW and S locations because these locations had higher 
values of Ts than the NE location. Since Ts gradients did exist between 
observation positions, sector horizontal homogeneity appeared to be 
influenced by the direction from which the air mass entered the area. 
The NW and NE positions, which experienced the largest Ts gradient, had 
more homogeneous ducting with northerly flow than westerly flow. The 
gradient between S and NE was small and resulted in consistently homo- 


geneous conditions in all three cases. 


9/7 





Case I illustrated the effect wind plays on Z* values. Both 
NE and NW had high relative humidity, but NW had a greater air-sea 
temperature difference, resulting in a larger moisture flux. However 
the stronger winds in the eastern region compensated for the lower 
moisture flux producing nearly similar duct heights. A comparison of 
the S and NE stations, which had similar winds and air-sea temperature 
differences, Showed Significantly different duct heights. The import- 
ant parameters in this case were relative humidity and sea surface 
temperature which produced larger surface moisture fluxes. 

Case II exhibits the effects due to incursions of cold, dry air 
with moderate to fresh winds. Satellite pictures (Figures 50, 51) 
indicated low level instability under windy and drying conditions. How- 
ever, Significant afternoon heating sharply reversed the air-sea tempera- 
ture difference causing the ducts to strengthen. Despite the afternoon 
variation, this case had the most consistent horizontally homogeneous 
ducting, occurring 44% of the time, as well as the most equal sector 
homogeneity of the three cases. The strongest ducts were toward the 
northwest where the air was the driest. The strong zonal flow in- 
hibited the driest from reaching the southern region where relative 
humidity varied little and air-sea temperature differences were small. 

Case III illustrates the effect of rapidly moving depressions 
on the formation and duration of ducts. The strongest areal ducts 
formed with the passage of a cold front but lasted only a short period 
before the next upstream weather system weakened the duct in the 


western region. The occurrence of ducts in the NE lagged the 
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occurrence of similar ducts at the NW position by six hours. 

A composite of the frequency of occurrence of areal and sector 
homogeneous ducts and maximum duration is depicted in Figure 18. It 
was concluded that areal homogeneity was a probable occurrence only 35% 
of the time for the three periods studied, with a maximum likely dura- 
tion of 12 hours. The NE sector appeared to maintain the best horizon- 
tal homogeneous conditions. The sector also appeared to exhibit good 
homogeneity but of short duration. The NW sector exhibited the least 
homogeneous conditions but when homogeneity occurred, it lasted up to 


24 hours. 


a2 


s 





81% (21) 


44% 75% 
(15) (30) 





a 54% (12) 





AVG. 65% (21} 





Figure 18. Percentage of areal and sector horizontal 
homogeneity of ducting and maximum duration for three 
cases and the average. 
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VII. SUMMARY AND CONCLUSIONS 


Unique temporal and spatial descriptions of evaporation duct 
features in the Northeast Atlantic were described from observational 
data taken from four ships spaced in a fixed array. Evaporation duct 
heights were determined from accurately measured bulk surface data 
utilizing Monin-Obukhov scaling in flux-profile relationships. Discern- 
ible synoptic features from surface pressure fields and DMSP satellite 
pictures were air mass histories, trajectories and horizontal variations 
(frontal zones). Synoptic scale spatial variations of the sea surface 
temperature were also available and were as important in the spatial 
variation as the atmosphere. 

Evaporation duct features established for this data set were: 

1. The mean Z* values were less than those reported in available 
climatologies based on weather ships’ observations (4.2m versus 
7.0m). 

2. Z* values at or above 8-10 meters occur frequently in conjunction 
with outbreaks of cold, dry air from the north or northwest and 
have durations of 9 to 12 hours. 

3. Air-sea temperature differences and relative humidity were the 
important surface layer parameters in establishing tactically 
Significant values and the former was the most important. 

4. If horizontal homogeneity were defined as representatively uniform 
Z* values over 200 km distances, such conditions occurred only 
35% of the times when significant duct heights existed, and for 
durations of up to 12 hours. 

Conclusions relative to synoptic scale features most closely asso- 
Ciated with these near-surface features are: 

1. The primary synoptic influence on the duct height is the air-mass 


trajectory (due to its role in determining the moisture content) 
and the air-surface temperature difference. North and northwest 
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trajectories have the highest associated Z* values. 

2. The primary feature influencing horizontal homogeneity over dis- 
tances of 200km is the sea surface temperature, and its influence 
depends on the air mass trajectory. 

5. Light wind conditions modify the above specified synoptic feature- 
surface duct feature relationships. The modifications occur over 
smaller temporal and spatial scales than those of the prevailing 
fie eine S|. 

A concluding assessment is that tactically significant features of 
the evaporation duct in the study area and for the study period can be 
established from synoptic scale descriptions. Further exploitation of 
available data should emphasize satellite imagery which could delineate 
surface temperature and air mass characteristics. 

A recommendation for further efforts is that all available surface 
data from the JASIN project be used in order to extend and substantiate 
feature descriptions and relationships identified in this study. Since 
the start of this study, additional surface layer data from buoys and 
other participating ships have been compiled by the organizing institu- 
tions and made available on computer tapes. These data would increase 


the observations used by two orders of magnitude, thereby improving the 


descriptions of both the temporal and spatial features. 
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APPENDIX A. 





TABLE A-I 
Julian (year) Day and Date 

Date Bay Day 
July 12 193 244 
13 194 245 
14 195 246 
1S 196 247 
16 197 248 
17 198 249 
18 199 250 
19 200 251 
20 20k 252 
1 202 253 
22 203 254 
23 204 255 
24 205 256 
25 206 257 
26 207 | 258 
28 209 259 
29 210 260 

30 211 

31 Pgs 
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Figure 19. (contd.) DAILY WEATHER MAPS 
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Figure 20, 500 mb (Top) and SFC Pressure Maps (Bottom), 
1200 GMT 19 Aug 78. 


69 








Pioume Zi, 500 
os Ate oa mb (Top) and SFC Pressure (Bottom) Maps, 1200 GMT, 


70 











500 mb (Top) and SFC Pressure (Bottom) Maps, 1200 GMT, 
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Figure 23. 500 mb (Top) and SFC Pressure (Bottom) Maps, 1200 GMT, 
25 Aug 78. 
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APPENDIX B 


Defense Meteorological Satellite Program (DMSP) imagery (visual 
and infrared) of the Northeast Atlantic and high-resolution imagery 


of the JASIN area corresponding to periods of significant ducting. 
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Figure 25. DMSP infrared imagery, 7 August, 10352 GMT. 
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Figure 28. DMSP visual imagery, 22 August, 1110 GMT. 





Figure 29. DMSP infrared imagery, 22 August, 1110 GMT. 
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Figure 31. DMSP infrared imagery, 23 August, 1221 GMT. 
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Figure 33. DMSP infrared imagery, 24 August, 1138 GMI. 
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Figure 34. DMSP visual imagery, 25 August, 1054 GMT. 
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Figure 35. DMSP infrared imagery, 25 August, 1054 GMT. 
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Figure 36. DMSP visual imagery, 
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Figure 37. DMSP infrared imagery, 28 August, 1235 GMI. 
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Figure 38. DMSP visual imagery, 29 August, 1150 GMT. 





Figure 39. DMSP infrared imagery, 29 August, 1150 GMI. 
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Figure 40. DMSP visual imagery, 30 August, 1107 GMT. 
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Figure 41. DMSP infrared imagery, 30 August, 1107 GMI. 
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DMSP infrared imagery, 31 August, 1218 GMT. 
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DMSP visual imagery, 31 August, 1218 GMT. 
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Figure 45. 





DMSP [Reamagery, 1 Sept, 1135 GMI. 
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imagery and surface data. 
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Figure 45. High-resolution DMSP IR imagery and surface data. 
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Figure 46. High-resolution DMSP visual imagery and surface data. 
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Figure 47. High-resolution DMSP IR imagery and surface data. 
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Figure 48. High-resolution DMSP visual imagery and surface data. 
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Figure 49. High-resolution DMSP IR imagery and surface data. 
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Figure 50. High-resolution DMSP visual imagery and surface data. 
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Figure 51. High-resolution DMSP IR imagery and surface data. 
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Figure 53. High-resolution DMSP IR imagery and surface data. 
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Fioure 54. High-resolution DMSP visual imagery and surface data. 
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Figure 56. High-resolution DMSP visual imagery and surface data. 
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Figure 57. High-resolution DMSP IR imagery and surface data. 
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Figure 59. High-resolution DMSP IR imagery and surface data. 
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